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O Abstract: We study several aspects of the kinetic approach to sterile neutrino production 

■ via active-sterile mixing. We obtain the neutrino propagator in the medium including 

self-energy corrections up to 0(G 2 F ), from which we extract the dispersion relations and o 
damping rates of the propagating modes. The dispersion relations are the usual ones 



in terms of the index of refraction in the medium, and the damping rates are T\(k) = 
Taa(fc) cos 2 m (k); ^(fc) = T aa (k) sin 2 9 m (k) where r aa (k) tx Gp kT 4 is the active neutrino 
scattering rate and 6 m {k) is the mixing angle in the medium. We provide a generalization 
of the transition probability in the medium from expectation values in the density matrix: 



and study the conditions for 



P a ^ s (t) = e -rit + e -r 2 t _ 2e -i(r 1+ r 2 )t cos (AEt) 

its quantum Zeno suppression directly in real time. We find the general conditions for 
quantum Zeno suppression, which for m s ~ keV sterile neutrinos with sin 28 < 10~ 3 may 
only be fulfilled near an MSW resonance. We discuss the implications for sterile neutrino 
production and argue that in the early Universe the wide separation of relaxation scales 
far away from MSW resonances suggests the breakdown of the current kinetic approach. 
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1. Introduction 



Sterile neutrinos, namely weak interaction singlets, are ubiquitous in extensions of the stan- 
dard model|, |, |, I and are emerging as plausible cold or warm dark matter candidates H 
H 0, H [| [H], [Tl], [l^, [l3|, [l4|, |l5|, 0, |l8], as potentially important ingredients in stellar col- 



lapse and supernovae[16, 19] and in primordial nucleosynthesis |20|, 21]. Sterile neutrinos 
with masses in the range ~ keV may also provide an explanation of pulsar "kicks" via 



asymmetric neutrino emission [22, 23]. 



The MiniBooNE collaboration [24] has recently reported results in contradiction with 
those from LSND[25, 2£] that suggested a sterile neutrino with Am 2 ~ 1 eV 2 scale. Al- 
though the MiniBooNE results hint at an excess of events below 475 MeV the analysis dis- 
tinctly excludes two neutrino appearance-only from z/„ — > v e oscillations with a mass scale 



Am 2 ~ 1 eV 2 , perhaps ruling out a light sterile neutrino. However, a recent analysis [27] 
suggests that while (3 + 1) schemes are strongly disfavoured, (3 + 2) neutrino schemes pro- 
vide a good fit to both the LSND and MiniBooNE data, including the excess of low energy 
events, because of the possibility of CP violation in these schemes, although significant 
tension remains between appearance and disappearance experiments. 



Sterile neutrinos as dark matter candidates would require masses in the keV range (5|, 
H @> H % 10) 11) 13 1 [H) 12) 18], and their radiative decay would contribute to the X-ray 
background H ^8|. Analysis from the X-ray background in clusters provide constraints on 
the masses and mixing angles of sterile neutrinos Jl3|, pS| , |3(i| , and recently it has been 
suggested that precision laboratory experiments on (3 decay in tritium may be sensitive to 
~ keV neutrinos 

Sterile neutrinos couple to standard model active neutrinos through an off diagonal 
mass matrix, therefore they are produced via active-sterile mixing. In the hot and dense 
environment of the early Universe when the scattering rate of active neutrinos off the 
thermal medium is large, namely a short mean free path, there is a competition between 
the oscillation length and the mean free path. It is expected that when the oscillation 
length is much larger than the mean free path, the active to sterile transition probability is 
hindered because rapid scattering events "freeze" the state to the active flavor state. This 
phenomenon receives the name of quantum Zeno effect or Turing's paradox, studied early 
in quantum optical coherence [33] but considered within the context of neutrino oscillations 
in a medium in references [34, 35, 36]. Pioneering work on the description of neutrino 
oscillations and decoherence in a medium was cast in terms of kinetic equations for a flavor 
"matrix of densities" [p^ l or in terms of 2 x 2 Bloch-type equations for flavor quantum 
mechanical states 0, j38| . A general field theoretical approach to neutrino mixing and 
kinetics was presented in [35, |36|, S, [l6| , however, while such approach in principle yields 
the time evolution of the distribution functions, sterile neutrino production in the early 
Universe is mostly studied in terms of simple phenomenological rate equations |||, 
39, 41, E3]. An early approach 



40 



relied on a Wigner-Weisskopf effective Hamiltonian for 
the quantum mechanical states in the medium, while numerical studies of sterile neutrinos 
as possible dark matter candidates (8|, [|!| rely on an approximate approach which inputs 
an effective production rate in terms of a time averaged transition probability p^, El] and 



relies on the following semiphenomenological rate equation [40, 41, 42, Ej^, || 

d 



dt 



f s (p,t) « T(a -> s;p) [f a {p;t) - f 8 (p;t)] 



fl.l) 



where f a)S are the distribution functions for active (a) and sterile (s) neutrinos, d/dt is the 
total time derivative including the redshift of momenta through the expansion in the early 
Universe and T(a — > s;p) is an effective reaction rate. It is determined to be|4(], ^l[ 

r„ 



Fa 



s;p) 



i(p)/p 
2 \ °" 



1.2 



where T aa (p) ~ G 2 F pT A is the active neutrino reaction rate and (P a ^ s ) is a time average 



of the active-sterile transition probability in the medium which in reference pl| is given by 

(p a _ s )=^y°°p a _ s (t)di (i.3) 

with P a ^ s the usual quantum mechanical transition probability but exponentially damped 



by a decoherence factor [41] 



Pa->s(t) 



sin 2 26L 
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[1 - cos(AM)] 



(1.4) 
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where AE, 8 m are the oscillation frequency and mixing angle in the medium respectively 
and Tdec = is the decoherence time scale. Hence the rate that enters in the kinetic 

equation (|1.1| ) is given by [41] 

F(a - s; p) = l J (1-5) 



1 + 



2ABfp)V 



Taa(p) ; 

where m (p), AE(p) are the mixing angle and active-sterile oscillation frequency in the 
medium respectively. The quantum Zeno paradox is manifest in the ratio 2AE(p) /T aa (p) in 
( |1.5| ): for a relaxation time shorter than the oscillation time scale, or mean free path smaller 
than the oscillation length, T aa (p) S> AE(p) and the active-sterile transition probability 
is suppressed, with a concomitant reduction of the sterile production rate in the kinetic 
equation (1.1). Most studies 42] of the production of sterile neutrinos via active-sterile 
mixing rely on the kinetic description afforded by equation ( |1.1| ). 

A field theoretic approach to sterile neutrino production near a MSW resonance which 
focuses primarily on the hadronic contribution and seemingly yields a different rate has 
been proposed in reference [44], and more recently it has been observed that quantum Zeno 
suppression may have important consequences in thermal leptogenesis]^]. 

Questions and goals: Recently we have studied the non-equilibrium aspects of oscil- 
lations and damping in a model of mesons that effectively describes the dynamics of mixed 
neutrinos in a medium in thermal equilibrium |46|]. In the case of two species of "neutrinos" 
this study reveals that there are two propagating modes in the medium, whose dispersion 
relations feature the index of refraction correction from forward scattering similar to those 
for neutrinos in the medium but also two different damping rates which are determined by 
the imaginary part of the self-energy correction evaluated on the mass shell. For the case 
of two mixed neutrinos in a medium it is natural to expect that the imaginary part of the 
self-energy corrections evaluated on the mass shell (dispersion relations) yield two different 
damping rates. Thus the results of ref.[^6| lead us to expect that for one active and one 
sterile mixed neutrinos the propagating modes in the medium feature two different damp- 
ing rates and this observation motivates the first question: How does the active-sterile 
transition probability P _> s (i) account for two damping scales?, namely why does the result 
( |1.4| ) feature only one damping scale?. This question leads to the related second ques- 
tion: How to generalize the concept of a transition probability to the case of propagation 
in a medium?. The usual transition probability is based on the evolution of single particle 
quantum mechanical wavefunctions which are linear superpositions of the eigenstates of 
the Hamiltonian. The statistical description of a medium does not rely on single particle 
wavefunctions or quantum mechanical states but on the quantum density matrix. Therefore 
the concept of the active-sterile transition probability P a -> s (t) in a medium must be gener- 
alized in terms of the quantum density matrix. As discussed above the final expression for 
the effective sterile production rate (F5) exhibits the Quantum Zeno suppression whenever 
F a a(p) 3> which has been argued to be the case at high temperature^. From the 

quantum field theory perspective this possibility is puzzling for the following reason: at 
high temperature the difference in the oscillation frequencies AE(p) is determined by the 



index of refraction correction from forward scattering [47] in the medium. This is deter- 



mined by a one- loop contribution [47] and is formally of order Gf, whereas the interaction 
rate r aa arises from an absorptive part of the self-energy and in the effective Fermi's field 
theory is at least of two- loop order, formally of order G F . Therefore from the field theo- 
retical perspective quantum Zeno suppression requires a competition of terms of different 
order in the perturbative expansion in Fermi's effective field theory [ 48 1. This observation 



brings us to the third question: considering and active neutrino with standard model 
interactions, can quantum Zeno suppression be manifest at high temperatures within the 
regime of validity of the perturbative expansion?. 

The emerging cosmological and astrophysical importance of sterile neutrinos motivates 
a deeper scrutiny of the current approach to the dynamical aspects of their production 
based on the rate equation ( |1.1| ) with the effective rate([1.5|). While at this stage we take 
this description based on ( |1.1[ - |1.3[ ) for granted, our goal is to address the three questions 
enunciated above within the quantum field theory of mixed neutrinos with standard model 
interactions in the medium, and in so doing we scrutinize the reliability of this approach. 
Our goals in this article are: i: to provide a quantum field theoretical understanding of 
the dispersion relations and damping rates of the two propagating modes (quasiparticles) 
in the medium, ii: to provide a generalization of the active-sterile transition probability 
in real time in the medium and a reassessment of the time averaged transition probability 
(P a _> s ) directly from the non-equilibrium time evolution of the full density matrix and iii: 
to scrutinize the possibility of quantum Zeno suppression within the realm of validity of 
perturbation theory with standard model interactions for the active neutrino. 

Main results: 



We consider one active and one sterile neutrino[21] to highlight the main conceptual 



aspects. Unlike most treatments in the literature that study the dynamics in terms of Bloch- 



type equations for a 2 x 2 flavor density matrix ]p4|, |4l| , 14], we study the full quantum field 
theoretical density matrix. A main advantage of studying the time evolution of the density 
matrix directly within the quantum field theory context is that we obtain the neutrino 
propagator which yields the quasiparticle dispersion relations and damping rates in the 
medium. Furthermore, the time evolution of the quantum field density matrix allows to 
study the non-equilibrium dynamics of neutrino mixing and propagation as an initial value 
problem from which we obtain the dispersion relation of the correct quasiparticle modes 
in the medium, their damping rates (widths) and the generalization of the active-sterile 
transition probability. These are all determined by the neutrino propagator in the medium 
which includes self-energy corrections up to two loops 0(G F ) in the standard model weak 
interactions. Our main results are the following: 

• There are two quasiparticle propagating modes, their dispersion relations are the 



usual ones in terms of the index of refraction in the medium|47] plus perturba- 
tive radiative corrections of (D(G F ) and their damping rates are given by Flip) = 
r aa (p) cos 2 6 m (p) ; r 2 (p) = T aa (p) sin 2 9 m (p) where T aa (p) oc CppT 4 is the active 
neutrino scattering rate in the absence of mixing, and 9 m (p) the mixing angle in the 
medium. We provide a physical interpretation of these different quasiparticle relax- 



ation rates and argue that these must naturally be correct in agreement with the 
fact that sterile neutrinos are much more weakly coupled to the plasma than active 
neutrinos in the regimes far away from an MSW resonance. 



We generalize the concept of the active-sterile transition probability P a ^, s (t) in the 
medium from expectation values of the active and sterile neutrino field operators 
in the full quantum density matrix. This is achieved by furnishing an initial value 
problem via linear response: the density matrix is initialized to feature an non- 
vanishing expectation value of the active neutrino field, but a vanishing expectation 
value of the sterile field. Upon time evolution a non-vanishing expectation value of 
the sterile neutrino field develops from which we extract unambiguously the transition 
probability. This formulation directly inputs the neutrino propagator with self-energy 
corrections up to 0{G 2 F ). We find 

Pa ^ si t) = [e-™ + e"™ - 2e-*( r i+r 2 )t cos (A£7i)J (1.6) 

This expression identifies the decoherence time scale for suppression of the interfer- 
ence term Tdec = 2/(^1 + ^2) = 2/r aa . Although the oscillatory interference term in 
the active-sterile transition probability is suppressed on the decoherence time scale 
Tdec, far away from an MSW resonance the relevant time scale for suppression of 
P a _> s (t) is determined by the smaller of the relaxation scales for the quasiparticles 
or Tzip). If the effective sterile production rate ( |1.2j ) is computed by inserting 
the result (O]) into the time average (1.3) the result for the effective production 



rate is enhanced far away from MSW resonances compared to that given by (|1.5|). 
However we argue that the widely different damping rates T12 suggest a breakdown 



of the simple rate equation (1.1) in these regimes in the early Universe. 



We provide a real time interpretation of the quantum Zeno suppression based on the 



generalization of the active-sterile transition probability (1.6). The complete general 
conditions for quantum Zeno suppression of the active-sterile transition probability 
are found to be: 

— a) the active neutrino scattering rate much larger than the oscillation frequency 
T aa (p) » AE(p), 

— b) The relaxation rates of the propagating modes must be approximately 
equal. In the case under consideration with Ti(p) = T aa (p) cos 2 6 m (p) ; ^(p) = 
r a a (p) sin 2 m (p) this condition determines an MSW resonance in the medium. 

Although these conditions are general, for sterile neutrinos with m s ~ keV and 
sin 20 < 1CU 3 and standard model interactions for the active neutrino, we find that 
they may only be fulfilled near an MSW resonance at Tmsw ~ 215 MeV, but a firm 
assessment of such possibility requires to include 0{G 2 F ) corrections to the index of 
refraction and a deeper assessment of the perturbative expansion. Far away from the 
resonance either at high or low temperature there is a wide separation between the two 



relaxation rates of the propagating modes in the medium. In these cases the transition 
probability reaches a TYidxiTRUTfi on the decohcrcficc time scale Td ec = 2/T aa , and is 
suppressed on a much longer time scale determined by the smaller of the damping 
rates. Even for T^AE^k) <C 1, which in the literature |34], |8[] is taken to indicate 
quantum Zeno suppression, we find that the transition probability is substantial on 
time scales much longer than Td ec if T± and I^ are widely separated. 

Section (|2|) provides a study of the time evolution of the full quantum field theory 
density matrix, and the equations of motion for expectation values of the neutrino fields. 
In this section we obtain the dispersion relations and widths of the propagating modes 
(quasiparticles) in the medium up to second order in the weak interactions. In section @ 
we introduced the generalized transition probability in the medium from the time evolution 
of expectation values of neutrino field operators in the density matrix. In this section we 
discuss in detail the conditions for the quantum Zeno effect, both in real time and in 
the time-averaged transition probability and the possibility of quantum Zeno suppression 
within the realm of validity of perturbation theory in standard model weak interactions. In 
section (|J) we discuss the implications of our results for the production of sterile neutrinos 
in the early Universe. In this section we argue that the in the early Universe far away from 
an MSW resonance the wide separation of the damping scales makes any definition of the 
time averaged transition probability ambiguous, and question the validity of the usual rate 
equation to describe sterile neutrino production in the early Universe far away from MSW 
resonances. Section (q) presents our conclusions . 



2. Quantum field theory treatment in the medium 
2.1 Non-equilibrium density matrix 

In a medium the relevant question is not that of the time evolution of a pure quantum 
state, but more generally that of a density matrix from which expectation values of suitable 
operators can be obtained. 

In order to provide a detailed understanding of the quantum Zeno effect, we need a 
reliable estimate of the dispersion relations and the damping rates of the propagating modes 
in the medium which are determined by the complex poles of the neutrino propagator in 
the medium. 

In this article we obtain these from the study of the real time evolution of the full 
density matrix by implementing the methods of quantum field theory in real time described 
in references p^, |50| , |5l| , 52, [5^, 46]. This is achieved by introducing external (Grassmann) 



sources that induce an expectation value for the neutrino fields. Upon switching off the 
sources these expectation values relax towards equilibrium and their time evolution reveals 
both the correct energy and the relaxation rates [f33|, [46]| . The main ingredient in this 
program is the active neutrino self-energy which we obtain up to second order in the 
standard model weak interactions. 

We consider a model of one active and one sterile Dirac neutrinos in which active- 
sterile mixing is included via an off diagonal Dirac mass matrix and the active neutrino 



only features standard model weak interactions. The relevant Lagrangian density is given 
by 



£ = £l + £ Ia , (2.1) 

where 

Cl = V(i$I-M)v, (2.2) 

with v being the neutrino doublet 



v a 



(2.3) 



and a, s refer to the flavor indexes of the active and sterile neutrinos respectively. 

The neutrino fields are four component Dirac spinors v = vr + vl with R, L the right 
and left handed components, both for a, s. The mass matrix in ( |2.2[ ) is of the Dirac type: o 
^Mi/£ + h.c. 

The Dirac mass matrix M is given by 



(2.4) 

o 
^1 



It can be diagonalized by the unitary transformation that takes flavor into mass eigenstates, 
namely 

(::H>(:)' <-> o 

with the unitary transformation given by the 2x2 matrix 

/ cos 9 rint\ ° 
\ — sin a cos t) I 

In this basis the mass matrix is diagonal 

M=( M ' ° V 

with the relation 

m aa = cos 2 9 Mi + sin 2 6M 2 ; m ss = sin 2 9 Mi + cos 2 9M 2 ; m as = ^(M 2 - Mi) sin 26 , 

(2.8) 

where 6 is the vacuum mixing angle. The Lagrangian density £j a describes the weak 
interactions of the active neutrino v a with hadrons or quarks and its associated charged 
lepton. Leptons, hadrons or quarks reach equilibrium in a thermal bath on time scales 
far shorter than those of neutrinos, therefore in what follows we assume these degrees of 
freedom to be in thermal equilibrium. Furthermore, in our analysis we will not include 



(2.7) 



1 



-7- 



the non-linearities associated with a neutrino background, such component requires a full 
non-equilibrium treatment and is not germane to the focus of this study. The Lagrangian 
density that includes both charged and neutral current interactions can be written in the 



form 1 36, 16 



£ia = ^ \O a Lv a + h. C .+V al ^Lv a ] (2.9) 

where L = (1— 7 5 )/2, the current includes both charge and neutral current contributions 
from the background in thermal equilibrium. For example the following charged and neutral 
current contributions to the effective Lagrangian (taking the active to be, for example, the 
electron neutrino): 



V2 



A[V e ^Le] [e^Lu e ] + 2V e ^[e lfl (g v - g Al *)e\Lv e } , (2.10) 



the first term (from charged currents) can be Fierz-rearranged to yield the form of the 
second term in fl2,9|) . This term yields a contribution to the index of refraction in the 
medium[[|7]], O a describes the charged current interaction with hadrons or quarks and the 
charged lepton, for example [36, ||] O a = j^L-ip^^^Cy — Ca7 5 )Vv I n the case of all 
active species the neutral current contribution to is the same for all flavors (when the 
neutrino background is neglected), hence it does not contribute to oscillations and the 
effective matter potential. In the case in which there are sterile neutrinos, which do not 
interact with the background directly, the neutral current contribution does contribute to 
the medium modifications of active-sterile mixing angles and oscillations frequencies. 

To study the dynamics in a medium we must consider the time evolution of the density 
matrix. While the usual approach truncates the full density matrix to a 2 x 2 "flavor" 
subspace thus neglecting all but the flavor degrees of freedom, and studies its time evolution 
in terms of Bloch-type equations [II]], our study relies instead on the time evolution of 
the full quantum field theoretical density matrix. 

The full density matrix describes a statistical ensemble of neutrinos and charged lep- 
tons, quarks or hadrons, these latter degrees of freedom are in thermal equilibrium and 
constitute the thermal bath. The fact that the density matrix describes charged leptons, 
quarks and or hadrons in statistical equilibrium will be used below (see eqns. 2.21, [2.22 ) 



when the correlation functions of these fields are obtained from ensemble averages in the 
density matrix. 

The time evolution of the quantum density matrix p is given by the quantum Liouville 
equation 

i^± = [H,p(t)] (2.11) 

where H is the full Hamiltonian with weak interactions. The solution is given by 

p(t) = e- iHt p{0)e iHt (2.12) 

from which the time evolution of observables associated with an operator A, namely its 
expectation value in the time evolved density matrix is given by 

(A(t)) =Ttp(t)A. (2.13) 



The density matrix elements in the field basis are given by 

p(^';t) = [ PW^|e-^ t |^)p(^^ / ;O)(0 / | e ^|^), (2-14) 



the matrix elements of the forward and backward time evolution operators can be handily 
written as path integrals and the resulting expression involves a path integral along a 
forward and backward contour in time. This is the Schwinger-Keldysh|4S, [3l|, 50, 52] 



formulation of non-equilibrium quantum statistical mechanics which yields the correct time 
evolution of quantum density matrices in field theory. Expectation values of operators are 
obtained as usual by coupling sources conjugate to these operators in the Lagrangian and 
taking variational derivatives with respect to these sources. This formulation of non- 
equilibrium quantum field theory yields all the correlation and Green's functions. Of 
primary focus is the neutrino retarded propagator 

S aP {x-x'-t-t') = -iQ{t-t') Trp{Q){^ a {x,t)^ p {x',t')}, (2.15) 

where the flavor indices a, (3 correspond to either active or sterile and ip a {x, t) is a neutrino 
field in the Heisenberg picture. The (complex) poles in complex frequency space of the 
spatio-temporal Fourier transform of the neutrino propagator yields the dispersion relations 
and damping rates of the quasiparticle states in the medium. It is not clear if this important 
information can be extracted from the truncated 2x2 density matrix in flavor space usually 
invoked in the literature and which forms the basis of the kinetic description ([O]), but 
certainly the full quantum field density matrix does have all the information on the correct 
dispersion relations and relaxation rates. 

A standard approach to obtain the propagation frequencies and damping rates of 
quasiparticle excitations in a medium is the method of linear response |5J]. An external 
source r\ is coupled to the field operator ip to induce an expectation value of this operator 
in the many body density matrix, the time evolution of this expectation value yields the 
quasiparticle dynamics, namely the propagation frequencies and damping rates. In linear 
response 

{tp a (x,t)) =Tr p(0)ip a (x,t) = - d 3 x'dt'S a/3 (x - x" ;t - t')n p {x\t') , (2.16) 



where S(x — x';t — t') is the retarded propagator or Green's function ( 2.15| ) and averages 



are in the full quantum density matrix. The quasiparticle dispersion relations and damping 
rates are obtained from the complex poles of the spatio-temporal Fourier transform of the 



retarded propagator in the complex frequency plane 1 54, pq |. For one active and one sterile 
neutrino there are two propagating modes in the medium. Up to one loop order 0(Gf) 
the index of refraction in the medium yields two different dispersion relations [[47]], hence we 
expect also that the damping rates for these two propagating modes which will be obtained 
up to 0(G 2 F ) will be different. This expectation will be confirmed below with the explicit 
computation of the propagator up to 0{G 2 F ). 

Linear response is the standard method to obtain the dispersion relations and damping 
rates of quasiparticle excitations in a plasma in finite temperature field theory[55|. The 



linear response relation ( |2.16| ) can be inverted to write 

S- 1 (iP(x,t)) = - V (x,t), (2.17) 

where 5 _1 is the (non-local) differential operator which is the inverse of the propagator, 
namely the effective Dirac operator in the medium that includes self-energy corrections. 
This allows to study the dynamics as an initial value problem and to recognize the quasipar- 
ticle dispersion relations and damping rates directly from the time evolution of expectation 
values of the field operators. This method has been applied to several different problems 
in quantum field theory out of equilibrium and the reader is referred to the literature for 



detailed discussions fl5q, 57, 5 



It is important to highlight that (ip(x,t)) is not a single particle wave function but 
an ensemble average of the quantum field operator in the non-equilibrium density matrix, 
namely an ensemble average. In contrast to this expectation value, a single particle wave 
function is defined as (l\ip(x,t)\0) where |0) is the vacuum and |1) a Fock state with one 
single particle. 

In the present case the initial value problem allows us also to study the time evolution 
of flavor off diagonal density matrix elements. Consider an external source r\ a that induces 
an initial expectation value only for the active neutrino field ip a , such an external source 
prepares the initial density matrix so that at t = the active neutrino field operator features 
a non-vanishing expectation value, while the sterile one has a vanishing expectation value. 
Upon time evolution the density matrix develops flavor off diagonal matrix elements and 
the sterile neutrino field ip s develops an expectation value. The solution of the equation of 
motion ( [2.17 ) as an initial value problem allows us to extract precisely the time evolution 



of (ip s ) from which we unambiguously extract the transition probability in the medium. 
2.2 Equations of motion in linear response 

The linear response approach to studying the non-equilibrium evolution relies on "adia- 
batically switching on" an external source r] that initializes the quantum density matrix 
to yield an expectation value for the neutrino field(s). Upon switching off the external 
source the expectation values of the neutrino fields relax to equilibrium. The real time 
evolution of the expectation values reveals the dispersion relations and damping rates of 
the propagating quasiparticle modes in the medium. These are determined by the poles of 
the retarded propagator in the complex frequency plane [54, 55]. 

The equation of motion for the expectation value of the flavor doublet is obtained by 



introducing external Grassmann- valued sources ??[pq, 57 



C s =Vri + fju, (2.18) 



shifting the field 



l/±=Va + ^ ; = ; (*a)=0, (2.19) 

for a = a, s, and imposing ) = order by order in the perturbation theory (5f], [57], pq| . 
Implementing this program we find the following equation of motion for the expectation 
value of the neutrino field ip a (x,t) = Tr[p(0)v a (x, t)] induced by the external field n 



(i 5 aP - M ap + Y%*£) ^(cc, t) + J d 3 x' J dt'^x-x', t - 1')^, t') = - Va (x, t), 

(2.20) 



The tadpole 



ytad = Gf^ o Tr ^ Q) j0(0) 0) (2 21) 

a/2 

describes the one-loop charged and neutral current contributions to the matter potential 
in the medium, and 

EZ\x - x<- t-t') = 1 -^- Trp(0) [o a (x, t)O a (x>, t') + O a (x>, t')O a [x, t)] . (2.22) 

The latter describes the two-loops diagrams with intermediate states of hadrons or quarks 
and the charged lepton, it is a fermionic correlation function in equilibrium and its spatio- 
temporal Fourier transform features an imaginary part that yields the relaxation rates of 



neutrinos in the medium. As shown in ref.[53], the spatial Fourier transform of the retarded 
self-energy can be written as 

Y, re \k, t-t')= - dk ImS(k, k ) e^*"*') . (2.23) 

7T J-oo 

The imaginary part ImE(fc,feo) evaluated on the mass shell of the propagating modes 
determines the relaxation rate of the neutrinos in the medium. Since only the active neu- 
trino interacts with the degrees of freedom in the medium, both self-energy contributions 
are of the form 



The initial value problem is set up as follows [53]. Consider an external Grassman 



valued source adiabatically switched on at t = — oo and off at t = 0, 

i] a (x, t) = rj a (x, 0) e et 9(-t) e -> 0+ . (2.25) 



It is straightforward to confirm that the solution of the equation of motion ( |2.20 ) for t < 
is given by 

il)p%t) =M^°) eet - ( 2 - 26 ) 

Inserting this solution for t < the equation of motion determines a relation between 
il)p{x., 0) and r] a (x, 0). For t > the equation of motion becomes an initial value problem 
with initial value given by ipg(x,0). For t > introducing spatial Fourier transforms and 
taking the Laplace transform, the equation of motion becomes (see ref. |5jJ for details) 



' • o 

Z7 s 



7 • kj 5 af3 - M a(3 + Y$f L + Z af3 (k, s) L\ ^(fc, a) = i ( 7 ° S a p + 0{G 2 F )) ^(k, 0) . 

(2.27) 



where ift, £ denote Laplace transforms with Laplace variable s. The Laplace transform of 
the retarded self energy admits a dispersive representation which follows from eqn.( 2~23| ), 
namely 1 53] 



E(A;,s) 



dko Im£(/c, ko) 



7T 



(2.28) 



Following ref. [|53| , we proceed to solve the equation of motion by Laplace transform as 
befits an initial value problem. 

In what follows we will ignore the perturbative corrections on the right hand side 
of (p^) since these only amount to a perturbative multiplicative renormalization of the 
amplitude, (see ref.|53|] for details). 



The chiral nature of the interaction constrains the self-energy to be of the form[47, 53] 



E tad L + s ^ s)L = { 7 u A(s; k) - 7 • kE(s; k) ) L 



(2.29) 



where the matrices A,B are of the form given in eqn. ( 2.24 ) with the only matrix elements 
being A aa ;B aa respectively. The dispersive form of the self-energy fl2,28| ) makes manifest 
that for s near the imaginary axis in the complex s-plane 



S(fe, s 



-iuj ± e) 



dko 



7T 



V 



Im£(fc, ko) 
ko — u) 



± iImS(A:, lo) 



(2.30) 



where V indicates the principal part. This result will be important below. 

The solution of the algebraic matrix equation ( p.27| ) is simplified by expanding the left 
and right handed components of the Dirac doublet ip in the helicity basis as 







v {h) (g, £(h) 





fe-E^^J ' * R= Ei r r 1 (2 ' 31) 

where the Weyl spinors are eigenstates of helicity a ■ k with eigenvalues h = ±1 and 
are flavor doublets with the upper component being the active and the lower the 
sterile neutrinos. 



Projecting the equation of motion ( 2.27 ) onto right and left handed components and 
onto helicity eigenstates, we find after straightforward algebra 



EC 

o 
o 

o 
o 
-J 

o 
o 
o 



[-{s 2 + k 2 )I + (is - hk){A{k; s) + hM{k\ a)) 



I 2 ] tp (h) (k,s) = i(is-hk)Iip {h) (k,0)-iM^ h) (jfc,0) 

(2.32) 



is + hk 
s 2 + k 2 



:^ h \k,s)+i^ h \k, 0) 



(2.33) 



where again we have neglected perturbatively small corrections on the right hand side of 
eqn. ( ggg ). 

It proves convenient to introduce the following definitions, 
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5M 2 = M\-Ml ; M 2 = -{Mf + 



Sh(k; s) = (is - hk)(A aa (k; s) + hB aa (k; s)) 
S h (k;s) 



A h (k;s) 



5M 2 



p h (k;s)= (cos 20 - A h (k;s)) 2 + sin 2 20 



cos 29$(k;s) 
29$(k;s) 



cos 29 - A h (k;s) 

Ph{k;s) 
sin 26 



sin 



Ph(k;s) 



(2.34) 
(2.35) 
(2.36) 

(2.37) 
(2.38) 

(2.39) 



in terms of which 



-(s 2 + k 2 )I + (is - hk)(A(k; s) +hM(k; s)) - M 2 

8M 2 f cos 29$ (k; s) -sm29$(k;s) 

— 2~PA( > s ) I _ sin2 ^)(fe ;5 ) -cos2^ ) (A;;s) 



- s 2 _fc 2 + -5 fc (fc,s)-M' 



(2.40) 



The solution of the equation (2.32) is given by 



^'(M) = sW(M) 



: C (h) (k, 0) + i(is - hk)I ip {h) (k, 0) 



(2.41) 



where the propagator §( h '(k, s) is given by 



S w (k, 



[a 2 (s,k)-(3 2 (s,k)] 



< ni.fl / / \{ cos29$(k;s) 



- sin 29$ (k; s) 
cos29$ (k;s) 



(2.42) 



o 
o 

o 
o 
-J 

o 
o 
o 



and we defined 



a h (k;s) 



(s 2 + k 2 ) + \s h (k;s)-M 2 



5M 2 

Ph(k;s) = — — ph(k;s) 



(2.43) 
(2.44) 



The real time evolution is obtained by inverse Laplace transform, 



ds 



si 



(2.45) 



where T is the Bromwich contour in the complex s plane running parallel to the imaginary 
axis to the right of all the singularities of the function (p(k, s) and closing on a large 
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semicircle to the left of the imaginary axis. The singularities of <p(k, s) are those of the 
propagator ( p. 42 ). If the particles are asymptotic states and do not decay these are isolated 



simple poles along the imaginary axis away from multiparticle cuts. However, in a medium 
or for decaying states, the isolated poles move into the continuum of the multiparticle cuts 
and off the imaginary axis. This is the general case of resonances which correspond to poles 
in the second or higher Riemann sheet and the propagator is a complex function with a 
branch cut along the imaginary axis in the complex s-plane as indicated by eqn. ( 2.3C| ). 



Its analytic continuation onto the physical sheet features the usual Breit-Wigner resonance 
form and a complex pole and the width determines the damping rate of quasiparticle 



excitations |5q, 57, 



It is important and relevant to highlight that the full width or damping rate is the sum 
of all the partial widths that contribute to the damping from different physical processes: 
decay if there are available decay channels, and in a medium the collisional width and or 
Landau damping also contribute to the imaginary part of the self-energy on the mass shell. 
The quasiparticle damping rate is one-half the relaxation rate in the Boltzmann equation 



for the distribution functions [p9|, 58] 



It is convenient to change the integration variable to s = —ioj + e with e — > + and to 
write the real time solution ( p. 45 ) as follows 



/OO J 
^ft h \k,s = -iu + e)e- i " t , (2.46) 
-oo 27r 

We focus on ultrarelativistic neutrinos which is the relevant case in the early Universe. 
Let us consider that initially there are no right handed neutrinos and only negative helicity 
are produced, namely h = — 1, and denoting the negative helicity doublet of expectation 
values as 

^•^fcgD' (2 - 47) 

where u a ^ s (k,t) now represent the expectation values of the negative helicity components 
of the neutrino fields in the density matrix. We find the expectation values at time t given 
by 

fv a {k,t)\ . f°° doj _ iuJt , , . (v a (%,G) 



i 



j_y^ +k)G{k;u) * ( , 48) 



where 



G(k;u) = S^(k,s = -iu + e) (2.49) 

and the integral in ( |2.48| ) is carried out in the complex lo plane closing along a semicircle 
at infinity in the lower half plane describing retarded propagation in time. 

In order to understand the nature of the singularities of the propagator, we must first 
address the structure of the self energy, in particular the imaginary part, which determines 



the relaxation rates. Again we focus on negative helicity neutrinos for simplicity. Upon 
the analytic continuation s = —icu + e for this case we define 



S(k, lj) = Sh=-i(k; s = -iuj + e) = (w + k)- Tr(7° - 7 • k)E oa (fc, s) 



s=—iui+e 



(2.50) 



From equation ( 2. 30]) which is a consequence of the dispersive form ( |2,28| ) of the self 
energy £ aa (A;,s) , it follows that 



S(k, to) = Sii(k, u) + iSj(k, uj) 



(2.51) 



where Srj are the real and imaginary parts respectively. The real part of the self energy 
determines the correction to the dispersion relations of the neutrino quasiparticle modes 
in the medium, namely the "index of refraction" , while the imaginary part determines the 
relaxation rate of these quasiparticles. 

2.3 The self-energy: quasiparticle dispersion relations and widths: 

Figure @) shows the one loop contributions of 0(Gf) including the neutral current tadpole 
diagrams which contribute to the in-medium "index of refraction" for one active species, 
and the two loop contribution of 0{Gp) with intermediate states of hadrons (or quarks) 
and the associated charged lepton, in the limit of Fermi's effective field theory. 
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Figure 1: Left: one loop contributions to the self energy, the diagrams in the second line yield 
S tad . These contributions are of 0(Gf). Right: two loops contribution of 0{G 2 F ) to the self-energy 
in Fermi's effective field theory limit, with internal lines corresponding to hadrons and the charged 
lepton, or alternatively quarks and the charged lepton above the QCD phase transition. 



o 



In a medium at temperature T the real part of the one-loop contributions to S(k, to) 
is of the form |7|, ||, H g 



S R (k,u) = (co + k)G F T 3 



L + 



T 



{aw + bk) 



(2.52) 



where L is a function of the asymmetries of the fermionic species and a, b simple coefficients, 



all of which may be read from the results in ref.[47, || 53]. Assuming that all asymmetries 
are of the same order as the baryon asymmetry in the early Universe L ~ 10~ 9 the term 



oc T/Myy in ( [2.52 ) for u> ~ k ~ T dominates over the asymmetry term for T > 3MeV[47, 
|53|| and in what follows we neglect the CP violating terms associated with the lepton 
asymmetry. 

The imaginary part to one loop order is obtained from a Cutkosky cut (discontinuity) 
of the diagrams with vector boson exchange shown on the left side in figure (^) and is 
determined by the processes W — ► lopa, Z — ► V a v a . Both of these contributions are expo- 
nentially suppressed at temperatures T <C Mw,z, hence the one- loop contributions to the 
imaginary part of S(k;uj) is negligible for temperatures well below the electroweak scale. 
The two loop contribution to the imaginary part is obtained from the discontinuity cut 
of the two loop diagram with internal hadron or quark and charged lepton lines in figure 
Some of the processes that contribute to the imaginary part in this order are for 
example neutron f3 decay n — > p + e + + V and its inverse, along with scattering processes 
in the medium. The imaginary part of the self-energy for these contributions on-shell is 
proportional to &; T 4 [[I?], 14] at temperatures T <C Mw- Therefore in this temperature 



range 

Si(k,ij) ~ (uj + k) G 2 F k T 4 . (2.53) 

The consistency and validity of perturbation theory and of Fermi's effective field theory 
for scales u,k,T <C M\y entail the following inequality 



Si(k,u) < S R (k,Lu) 



(2.54) 



For example near the neutrino mass shell for ultrarelativistic neutrinos with to ~ k, assum- 
ing L ~ 10~ 9 and discarding this CP violating contribution for T > 3MeV because it is 
subleading, we find 

Si(k,uj) 2 

Q-7Z — V ~ 9 ( 2 - 55 ) 
S R {k,u) 

with g the standard model weak coupling. This discussion is relevant for the detailed 
understanding when the quantum Zeno effect is operative (see section (|j) below). 
The propagator G(k;ui) for negative helicity neutrinos is found to be given by 



cow = £ 



1 



1 



a + (3 cos 29 m —fi sin 29 m 
—(3 sin 29 m a — (3 cos 29 r , 



a — (3 a + 

where we have suppressed the arguments for economy of notation, and defined 



(2.56) 



a = uj 2 -k 2 -M 2 + -S R (k,uj) + -S I (k,LO 



(3 



5M 2 



cos 29 



2 

S R (k,uj) 



5M 2 

°m = 9 { - l) {k,s = -iuj + e). 
The inequality ( p. 54 ) licenses us to write f3 consistently up to 0{G F ) as 



(2.57) 

(2.58) 
(2.59) 




vww 



Figure 2: Contributions to the imaginary part of the self energy, the vertical line represents a 
Cutkosky cut. Left: discontinuity from the one loop contributions to the self energy of 0(Gf), 
from the decay of vector bosons for example W — > W, Right: discontinuity from the two loops 
contribution of 0{G 2 F ), arising for example from n — > p + e + + V e or similar processes at the quark 
level. 

O 
-~J 



where 



p(k,u>) 



5M< 



p(k,u) — -Sj(k, uj) cos 29 rt 



O 
O 

o 



COS °~ SAP ) +Sm ° 



(2.60) 



(2.61) 



Equation ( 2.56| ) makes manifest that G[k\uj) is strongly peaked at the values of uj 
for which a = ±/3. These determine the position of the complex poles in the analytic 
continuation. In the relativistic approximation k 3> M\ 2 we find: 



o 
o 

N) 
o 
o 
-j 

o 
o 
o 



For a = (3: 



with 



Ei{k) 
T x (k) 



uji{k) = E x {k) - i 



ri(fc) 



2fc 



"TT2 ^ 

M + 



2 PV ' 2 



cos 2 6> m (fc) 



(2.62) 

(2.63) 
(2.64) 



For a 



-0: 



io 2 (k) = E 2 (k) - i 



r 2 (*0 



(2.65) 
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with 



E 2 (k) 

r 2 (fc) 



k + 
r 

- 1 - an. 



1 

(*) 



M' 



■ sin 



,5M 2 S fl (A;) 
-z~PW =— 



(2.66) 
(2.67) 



where 



p{k) 



(2.68) 



and T aa (k) is the standard model result for the scattering rate of the active neutrino 
speciesg, 0, 0, H 



Taa(fe) Sl(k,U = k) 1 



^Tr( 7 u - 7 • k)ImS aa (A ; , u; = k)~G 2 F T 4 k 



(2.69) 



2 2fc 4 

and B m (k) = #m~ ^(^j 5 = — ifc) is the mixing angle in the medium for negative helicity 
neutrinos of energy u ~ k in the relativistic limit. The relations ( [2.64 , Oi7| ) are the same 



as those recently found in reference [46]. 
Combining all the results we find 



v a (k,t) 
v s (k, t) 



e - iEl { k )t e - T -±^t I J l + cos28 m {k) -sm29 m (k) ^ + 



sin2# m (£;) 1 — cos 29 m (k) 



This expression can be written in the following more illuminating manner, 



e -iE 2 {k)t e -^M t 1 / 1 - cos 26 m {k) sin 29 m (k) 

2 I sin20 m (£;) l + cos26» m (fc) 



(2.70) 



o 
o 

N) 
o 
o 
-j 

o 
o 
o 



\Vs(k,t) J \ e -iE 2 (k)t e -^-t J yv s (k,0) J 

(2.71) 

where U(0 m (k)) is the mixing matrix ( |2.6D but in terms of the mixing angle in the medium. 

In obtaining the above expressions we have neglected perturbative corrections from 
wave function renormalization and replaced to + k ~ 2k thus neglecting terms that are 
subleading in the relativistic limit, and the imaginary part in u>, which although it is of 



0(Gp), yields the effective Wigner-Weisskopf approximation [46| 



2.4 Physical interpretation: 

The above results have the following clear physical interpretation. The active (a) and 
sterile (s) neutrino fields in the medium are linear combinations of the fields associated with 
the 1,2 quasiparticle modes with dispersion relations E\2(k) and damping rates Ti^k) 
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respectively, on the mass shell of the quasiparticle modes the relation between them is the 
usual one for neutrinos propagating in a medium with an index of refraction, namely 



v a (k,t) = cos9 m (k) vi(k,t) +sin9 m (k) 1*2$, t) (2.72) 
v s (k,t) = cos9 m (k) U2(k,t) — sm9 m (k) vi(k,t) . (2.73) 

These relations between the expectation values of flavor fields and the fields associated 
with the propagating quasiparticle modes in the medium are obtained from the diagonal- 
ization of the neutrino propagator on the mass shell of the quasiparticle modes. These are 
recognized as the usual relations between flavor and "mass" fields in a medium with an 
index of refraction. 

At temperatures much higher than that at which a resonance occurs (and for k ~ T) 
9 m (k) ~ 7r/2 then v a ~ 1/2, and the active neutrino features a damping rate T2 ~ r aa 
while the sterile neutrino v s ~ v\ with a damping rate T\ = T aa cos 2 9 m (k) <C r aa . In 
the opposite limit for temperatures much lower than that of the resonance and for very 
small vacuum mixing angles v a ~ v\ and features a damping rate T\ ~ r aa while v s ~ Vi 
with a damping rate T2 r^Q, sm 2 9 <§C Taa* Thus it is clear that m both limits the 
active neutrino has the larger damping rate and the sterile one the smallest one. This 
physical interpretation confirms that there must be two widely different time scales for 
relaxation in the high and low temperature limits, the longest time scale or alternatively 
the smallest damping rate always corresponds to the sterile neutrino. This is obviously 
in agreement with the expectation that sterile neutrinos are much more weakly coupled 
to the plasma than the active neutrinos for sm29 m (k) ~ 0. This analysis highlights that 
two time scales must be expected on physical grounds, not just one, the decoherence time 
scale, which only determines the suppression of the overlap between the propagating states 
in the mixed neutrino state. 

The solution (^?70| ) can also be obtained from an effective Schroedinger-like equation 



but with an effective Weisskopf-Wigner Hamiltonian[46], not for the quantum states, which 
are not decaying, but for the ensemble averages, 

-H-f^!!) (2-74) 
v s {k,t) J \v s {k,t) J 

where the Weisskopf-Wigner effective Hamiltonian is 
M\ 5M 2 ( cos2# -sin2fl\ 1 





«~ " ( k + -2k) 1 + -UT { - sin 20 -COS2.J-2* \ S "^ +iS '^>l 1^1- ("5) 

Although this form of the time evolution looks similar to the usual quantum mechanical 
one, we emphasize that this equation of motion is for the ensemble averages of the neutrino 
operators and Srj are obtained from the self-energy on the mass shell of ultrarelativistic 
neutrinos. This self energy is an ensemble average of correlation functions of the operators 
that describe the degrees of freedom in the medium. This time evolution is an a posteriori 
consequence of the field-theoretical analysis, which unambiguously reveals the time evo- 
lution. While the description afforded by the Schroedinger-like equation for the ensemble 



averages with a non-hermitian Weisskopf-Wigner Hamiltonian is simpler, it can only be 
rigorously justified through the detailed study presented above. 

The emergence of two time scales can also be gleaned in the pioneering work on ster- 



ile neutrino production of ref.[40| (see eqn. (10) in that reference) which in this refer- 
ence where obtained within a phenomenological Wigner-Weisskopf approximation akin to 
( 2.74 , 2.75 ). Our quantum field theory study based on the full density matrix and the 
neutrino propagator in the medium provides a consistent and systematic treatment of 
propagation in the medium that displays both time scales and provides a derivation of the 
effective Schroedinger-like evolution with the Weisskopf-Wigner Hamiltonian 46]. 



3. Quantum Zeno effect 

3.1 Real time interpretation and general conditions 

Consider a density matrix in which the expectation value of the active neutrino field is non- 
vanishing, but that of the sterile neutrino field vanishes at the initial time t = 0, namely 
v a (k,0) ^ ; f s (k, 0) = 0. Then it is clear from equation fl2.70| ) that flavor off-diagonal 
density matrix elements develop in time signaling that sterile neutrinos are produced via 
active-sterile mixing with amplitude 



u s (k,t) 



-sm26 m (k) 



e -iE-,{k)t e 



e -iE 2 (k)t t 



v a (k,0) (3.1) 



From the solution ( |3.1| ) we introduce the generalized transition probability in the medium 
from the expectation values of the neutrino fields in the density matrix, these are the 
transition probabilities between ensemble averages of one-particle states of the neutrino 
fields, 



v s (k,t) 



f«(M) 



sin 2 29 m (k) 



-r(k)t 



cosh(j(k)t) - cos(AE(k)t) 



(3.2) 



EE 

o 
o 

o 
o 
-J 

o 
o 
o 



where 



T(k) = l -(r 1 (k) + r 2 (k))=^l 

7 (fc) = i(ri(AO-r 2 (*o) = ^^cos2Mfc) 



AE(k) = E 1 (k)-E 2 {k) 



5M Z 
~2k 



p(k) 



For the analysis that follows it is more convenient to write ( ^.2[ ) in the form 



; (k;t) 



sin 2 29 m (k) 



+ e -r 2 (k)t _2e- r( - k)t cos(AE(k)t) 



(3.3) 
(3.4) 
(3.5) 



(3.6) 
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The first two terms are obviously the probabilities for the quasiparticle modes 1,2, while 
the oscillatory term is the usual interference between these but now damped by the fac- 
tor e _r ( fc )*. This form of the transition probability is remarkably similar to the familiar 
transition probability in Kq — Kq or Bq — Bq oscillations |30], 61]. However we emphasize 
that ( [3.6] ) is the generalized transition probability extracted from expectation values in the 
density matrix. 

We highlight that the decoherence time scale is precisely T~ 1 (k) = 2/T aa (k) as antic- 



ipated in references [34, 41], since the interference between the two quasiparticle modes is 
suppressed on this time scale. However, the total transition probability is suppressed on 
this time scale only if Y\ = T2 = T, namely near a resonance. In this case 



.(*) 



sim 20 m e 2 1 sm z 



AE 



(3.7) 



which is the result quoted in reference |41| (see eqn. (|1.4| )). Under these conditions quantum 
Zeno suppression occurs when T(k) 3> AE{k) in which case the decoherence time scale is 
much smaller than the oscillation time scale and the transition probability is suppressed 
before a —* s oscillations take place. 

However, far away on either side of the resonance, although the oscillatory interference 
term is suppressed on the decoherence time scale T -1 , the transition probability is not 
suppressed on this scale but on a much longer time scale, determined by the smaller of 
Ti t 2. Only when T± = T2 = T, namely 7 = both the coherence (oscillatory interference 
term) and the transition probability are suppressed on the decoherence time scale. This 
phenomenon is displayed in figures (HJ4| ) which show the transition probability as a function 
of time without the prefactor sin 2 26 m (k) for several values of the ratios 7/T; AE/T. 
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o 
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Figure 3: The transition probability P a _» s (i) (without the prefactor sin 2 2# m /2) vs. Ft. The 
figure depicts the cases cos29 m (k) — 0.98 and cos2# TO (fc) = respectively, both with AE/T = 
S M 2 p(k) / 2kT = 5. The scale for suppression of the oscillatory interference is 1/r in both cases. 



Even for T(k) S> AE(k), claimed in the literature |}4|, ^] to be the condition for 
quantum Zeno suppression, the transition probability is substantial on time scales much 
longer than T _1 if T\ and are widely separated, namely if \^f/T\ ~ 1. This situation is 
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depicted in figure @ . From this analysis we conclude that the conditions for quantum Zeno 
suppression of P a ^ s (t) are: i): F(k) 3> AE(k) and ii):) j(k) ~ 0, namely T\(k) ~ T2(k). 
These conditions are obtained directly from the time dependence of P a ^ s (t) without taking 
any time average. We then emphasize that it is not necessary to average the probability 
over time to recognize the criteria for the quantum Zeno effect, these can be directly gleaned 
from the time evolution of the probability as originally proposed 




Figure 4: The transition probability P a _ >s (t) (without the prefactor sin 2 2# m /2) vs. Tt. The left 
panel is for 7/r = 0.5, AE/T = 4, 1,0.1, the right panel is for 7/T = 0.2, AE/T = 4, 1,0.2. 



From the arguments in reference M], the effective sterile neutrino production rate is 
obtained from the average of the transition probability on the decoherence time scale T^ ec . 
Using the result fl3.6|) we find instead 



Pa- 



Pa^s(t)dt 



sin 2 26V 



(3.8) 



where r, 7 are given by eqns. (3.3,3.4) respectively. 

This expression features two remarkable differences with the result (|l.5| ) f4l|l : the extra 
terms (7/r) 2 in the numerator and 1 — (7/r) 2 in the denominator, both are consequence 
of the fact that the relaxation is determined by two time scales Ti,r2. Only when these 
scales are equal, namely when 7 = the result ( |1.5| ) often used in the literature is recovered. 
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y/r = 0.98 
AE/r =0.1 





y/r=o 
AE/T= o.i 



rt 

Figure 5: The transition probability P a ^ s (t) (without the prefactor sin 2 2# m /2) vs. Ft in the 
quantum Zeno limit r > A£ for the cases cos26» m (fc) = 0.98,0 and AE/T = SM 2 p(k)/2kT = 0.1. 
The ramp-up time scale is ~ 1 /Y\ ~ 1 /V. In the left figure the damping time scale is ~ 1 /T% ~ 50/T. 
The right figure displays the resonant case for which the damping and coherence time scale coincide, 
when the conditions for quantum Zeno suppression arc fulfilled. 



We note that the result ( |3.8[) is not singular since 1 — (7 2 /T 2 ) = sin 2 26 m , (see also the 
discussion in section (|4.1D ), 

This analysis leads us to state that the complete conditions for quantum Zeno suppres- 
sion of the transition probability are that both 7/r < 1 and AE/T <C 1. That these are 
indeed the correct necessary conditions for quantum Zeno suppression can be gleaned from 
figures (|j, ||) which display the transition probability (without the prefactor sin 2 29 m /2) as 
a function of time for several values of the ratios 7/r, AE/T without performing the time 
average. 



o 
o 

ISO 
O 
O 
-J 

o 
o 
o 



3.2 High and low temperature limits: assessment of the quantum Zeno condi- 
tion 

In order to establish when the quantum Zeno condition AE(k)/T aa (k) <C 1 is fulfilled we 
focus on the cases far away from resonances and, according to the exhaustive analysis of 
ref.p, H and the constraints from the X-ray background in clusters [29, 30], in the region of 
parameter space IkeV < m s < lOkeV , 1CT 10 < sin 2 29 < 10~ 6 . We consider T > 3MeV 
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for which we can neglect the CP violating asymmetry contribution in ( 2.52 ) assuming 
that it is of the same order as the baryon asymmetry L ~ 1CP 9 [47, 53]. In this regime 
5M 2 ~ m 2 , and from ( |2.52| ) we find 



5M 2 \ MeV 



???,.„. 



Taking k ~ T and m s ~ IkeV the MSW resonance Sji(k, k)/5M 2 = 1 occurs at 
Tmsw ~ 215 MeV (a more precise estimate yields T ~ 180MeV|, |). For T > T M sw 
corresponding to Sn(k, k)/5M 2 3> 1 the active sterile oscillation frequency becomes 



5M 2 /M S R (k,k) G F T^k 

L w 

From the result ( |2.69| ) for r aa (£;) we find in the high temperature limit T 3> Tmsw 



AE(k) = ~ ~ -Jg- (3.10) 



2AE(k)_ G F T'k ... 



T aa (k) G 2 F T 4 k Myy g 2 

where g is the weak coupling. We note that in the high temperature limit the ratio 
AE(k)/Y aa (k) becomes independent of T, k. This result is in agreement with the con- 
clusions in ref.Q. 

In the low temperature limit 3 MeV < T -C Tmsw it follows that Sji(k,k)/5M 2 <C 1 
and the active-sterile oscillation frequency is 

2k 

hence the ratio 



AE(k) ~ ^ (3.12) 



r^)~GP^~ 10 ^kevj ^MevJ J ^ 

which for A; ~ T can be simplified to 

At the MSW resonance T = Tmsw, cos2# m ~ 0, AE(k) = m 2 sin26 /2k and the ratio 
becomes 

2AE(k) 9 , 
- ~ 10 2 sin 29 < 1 3.15 

for 10~ 5 < sin 26 < 10 3 . Therefore at the MSW resonance cos29 m (k) ~ and both 
conditions for quantum Zeno suppression, j/Y <C 1 , AE/Y <C 1 are fulfilled. However, 
we point out that near the resonance the oscillation frequency becomes AE 
onZy if the second order corrections to the dispersion relations (real part of the poles) are 
neglected, therefore if Y S> AE a reassessment of the perturbative expansion including 



second order corrections to the dispersion relations is required, since F oc G 2 F . Therefore 
this analysis leads us to unambiguously conclude that with standard model interactions for 
the active neutrino quantum Zeno suppression is not realizable either at high temperatures, 
when the matter potential dominates or at very low temperatures where the mixing angle is 
close to the vacuum value. Such possibility may only emerge very near an MSW resonance, 
however for small mixing angle this case requires a thorough reassessment of the dispersion 
relations in the medium including corrections of 0{G 2 F ) to the oscillation frequency. 

3.3 Validity of the perturbative expansion: 

The quantum Zeno condition T aa {k) S> AE(k) requires a consistent assessment of the 
validity of the perturbative expansion in the standard model and or Fermi's effective field 
theory. The active neutrino scattering rate r aa oc G 2 F kT A is a two loops result, while to 
leading order in weak interactions, the index of refraction contribution to the dispersion 
relation Sji(k,u)) is of one-loop order |47], |53|| . In the high temperature limit when Sr » 
SM 2 ~ m 2 s the active-sterile oscillation frequency is 

AE(k) ~ (3.16) 

combining this result with equation ( |2.69| ) at high temperature or density where the index 
of refraction dominates over 5M 2 , it follows that 

A£(fc) \S R (k,k)\ 

T aa (k) ~ Si(k, k) 1 ' ; 

for k ~ T the perturbative relation ( |2,55[ ) states that this ratio is > l/<7 2 3> 1 where g is 
the weak gauge coupling. An opposite ratio, namely AE(k)/T aa <C 1 would entail that the 
two-loop contribution (T aa ) is larger than the one- loop contribution that yields the index 
of refraction Sr. Thus quantum Zeno suppression at high temperature when the index of 
refraction dominates the oscillation frequency necessarily implies a breakdown of the strict 
perturbative expansion. Such potential breakdown of perturbation theory in the standard 
model or Fermi's effective field theory in the quantum Zeno limit has been already observed 
in a different context by these authors in ref.|48|], and deserves deeper scrutiny. We are 
currently exploring extensions beyond the standard model in which neutrinos couple to 
scalar fields motivated by Majoron models, in these extensions the coupling to the scalar 
(Majoron) provides a different scale that permits to circumvent this potential caveat. We 
expect to report on our results in a forthcoming article|p2|. 



4. Implications for sterile neutrino production in the early Universe: 

4.1 Time averaged transition probability, production rate and shortcomings of 
the rate equation 



The effective sterile production rate proposed in ref.[41| and given by eqn. ( |l.2[) requires the 



average of the transition probability P a _> s (i) over the decoherence time scale. Hence, com- 



bining ( |1.3| ) with the transition probability in the medium ( ^3.2|j3T3 ,3.4) yields the following 
time averaged transition probability (compare to eqn. (| 



( 2AE(k) \ 2 



sin 

r a— >s 



sin 2 29 m {k) 



l + (^f?) 2 



(4.1) 



We have purposely kept the sin 2 29 m (k) in the numerator and denominator to highlight 
the cancelation between this factor arising from the transition probability in the numerator 
with the factor 1 — (7/r) 2 arising from the total integrated probability in the denominator. 
The factor cos 2 29 m {k) in the numerator and the sin 2 29 m (k) in the denominator are hall- 
marks of the presence of the two different relaxation rates Ti(k),T2(k), and are responsible 
for the difference with the result (1.5). The extra factor sin 2 2# m (£?) in the denominator 
signals an enhancement when 9 m (k) = 0, tt/2. In the case 9 m (k) ~ the relaxation rate 
^2{k) <C ^i(k) whereas for 9 m (k) ~ tt/2 the opposite holds, ri(fc) <C ^(fc). In either case 
there is a wide separation between the relaxation rates of the propagating modes in the 
medium and the longest time scale for relaxation dominates the time integral in ( [4.1[ ) . This 
is depicted in fig. (]5|). 



(4.2) 



The result of references ][il|, 42, |8|] is retrieved only near an MSW resonance for which 
cos2# m (fc) ~ 0, in this case the relaxation rates become the same and 7 = 0. However, 
accounting for both relaxation rates I\;r2 yields the new result ( f4.1|]4.2| ) which is generally 
very different from the usual one (|1.5|). 

The result ( |4.2| ) is in clear contradiction with the analysis in section ([T4|) wherein the 
physical interpretation of the damping rates identify the sterile degrees of freedom as very 
weakly coupled to the plasma both at high temperature {9 m ~ tt/2) and low temperature 
(9 m ~ 0), therefore should feature small production rates. Contrary to this expectation, 
taking the limit of sin2# m (£;) ~ in ( [4 . 2[) , still yields a non- vanishing sterile neutrino 



This is an important difference with the result in [41] wherein it was assumed that 
T\ = T2, in which case 7 = 0. For 9 m (k) ~ 0, tt/2, the ratio 7/r ~ 1 leads to an 
enhancement of the time averaged transition probability. The interpretation of this result 
should be clear. The probability P a ^ s (t) has two distinct contributions, the interference 
oscillatory term, and the non-oscillatory terms. When one of these non-oscillatory terms 
features a much longer relaxation time scale, it dominates the integrand at long time after 
the interference term has become negligible, as shown in figure (|5]). Therefore the time 
integral receives the largest contribution from the term with the smallest relaxation rate, 
this is the origin of the factor 1 — (7/r) 2 = sin 2 29 m {k) in the denominator, o 

Taking the kinetic equation that describes sterile neutrino production ( |1.1| ) along with 
the effective production rate ( |1.2[ ) at face value, the new result (|4.1|) for the average 
transition probability yields the effective production rate 



^1 



production rate despite the fact sterile neutrinos decouple from the plasma in this limit. 
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The origin of this puzzling result is the time averaged probability (P a ^ s ) QL3|) and not 

any ambiguity in the calculation of the relaxation rates or in the time dependence of the 



transition probability P a _, s {t). The time integral in the averaged expression (1.3) introduces 
a denominator sin 2 29 m (k) from the longest time scale, and it is this denominator that is 
responsible for the enhancement. Thus the unreliability of the result ( ^4.2[ ) is a direct 
consequence of using the time-averaged transition probability ( |1.3| ) in the rate equation 

The real time analysis presented above clearly suggests that far away from an MSW 
resonance when T\ and T2 differ widely, T _1 is not the relevant time scale for suppression 
of the transition probability, but the longest of and T^ 1 therefore the time averaged 



transition probability (1.3) cannot be the correct ingredient in the rate equation. A more 
suitable definition of the average transition probability under these circumstances should 
be 

poo 

{P a ^ s ) = T sm / P a ^ s (t)dt (4.3) 



J 

where T sm is the smallest of T\ t 2- In a non-expanding cosmology this would indeed be the 
correct definition of an average transition probability, however in the early Universe as the 
temperature diminishes upon cosmological expansion, T sm changes with time crossing from 
T\ over to T2 at the resonance and the alternative definition (|4.3j ) would imply a "rate" 
with a sliding averaging time scale that changes rapidly near an MSW resonance. One can 
instead provide yet another suitable definition of an averaged transition rate 



p -p poo 

<P a _ s > = — / Pa _> s (t) dt . 

1 1 + 1 2 Jo 



(4.4) 



When the two rates differ widely the prefactor always approximates the smaller one. Since 

Ti T2 T aa 



T l + T 2 4 



T 2 (k) 



(4.5) 



this definition would cancel the enhancement from the sin 2 29 m (k) in the denominator in 
( |4.1| ) (still leaving the cos 2 29 m (k) in the numerator), but it misses the correct definition 
of the average rate by a factor 2, namely by 100%, in the region of the resonance where 
Ti = r 2 = T aa /2. 

Obviously the ambiguity in properly defining a time averaged transition probability 
stems from the wide separation of the time scales associated with the damping of the 
quasiparticle modes, far away from an MSW resonance. Near the resonance both time 
scales become comparable and there is no ambiguity in the averaging scale. Complicating 
this issue further is the fact that in the early Universe these time scales are themselves time 
dependent as a consequence of the cosmological expansion and feature a rapid crossover 
behavior at an MSW resonance. 



4.2 Caveats of the kinetic description. 

It is important to highlight the main three different aspects at the origin of the enhanced 
production rate given by equation fl4.2|) in the high temperature regime, for 9 m (k) ~ tt/2: 



i) the assumption of the validity of the usual rate equation in terms of a time-averaged 
transition probability wherein the relevant time scale for averaging is the decoherence 
time scale 1/T, ii) the result of a complete self-energy calculation that yields two time 
scales which are widely different far away from an MSW resonance, in particular at very 
high and very low temperatures, iii) the generalized transition probability in the medium 
Q3.6| ) obtained from expectation values in the full density matrix, rather than the usual 
quantum mechanical expression in terms of single particle states. The real time study of 
the transition probability shows that the oscillatory interference term is suppressed on 
the decoherence time scale 1/T, but also that this is not the relevant time scale for the 
suppression of the transition probability far away from an MSW resonance. The transition 
probability actually grows during 1 /T reaches its maximum on this time scale and remains 
near this value for a long time interval between 1/T and 1/T sm where T sm is the smaller 
of ri,T2. The enhanced rate emerges when taking for granted the definition of the time- 
averaged probability in terms of the decoherence scale but including in this expression the 
correct form of the transition probability Q3.6| ). As discussed above, alternative definitions 
of a time-averaged rate could be given, but all of them have caveats when applied to sterile 
neutrino production in the early Universe. 

However, we emphasize, that the underlying physical reason for the enhancement does 
not call for a simple redefinition of the rate but for a full reassessment of the kinetic 
equation of sterile neutrino production. The important fact is that the wide separation of 
scales prevent a consistent description in terms of a simple rate in the kinetic equation, 
a rate implies only one relevant time scale for the build-up or relaxation of population, 
whereas our analysis reveals two widely different scales that are of the same order only 
near an MSW resonance. 

Kinetic rate equations are generally a Markovian limit of more complicated equations 
in which the transition probability in general features a non-linear time dependence . Only 
when the non-linear aspects of the time dependence of the transition probability are tran- 
sients that disappear faster than the scale of build-up or relaxation an average transition 
probability per unit time, namely a rate, can be defined and the memory aspects associ- 
ated with the time evolution of the transition probability can be neglected. This is not 
the case if there is a wide separation of scales, and under these circumstances the assump- 
tions leading to the kinetic equation ( |1.1| ) must be revised and its validity questioned, very 
likely requiring a reassessment of the kinetic description. This situation becomes even more 
pressing in the early Universe. In the derivation of the average probability in ref . [ 41 1 the 
rate T aa (denoted by tq in that reference) is taken as a constant in the time integral in the 
average. This is a suitable approximation if the integrand falls off in the time scale 1/T aa , 
since this time scale is shorter than the Hubble expansion time scale for T > 1 MeV. How- 
ever, if there is a much longer time scale, when one of the relaxation rates is very small, 
as is the case depicted in fig.(|5|), then this approximation cannot be justified and a full 
time-dependent kinetic description beyond a simple rate equation must be sought. 

Thus we are led to conclude that the simple rate equation (|l]l]) based on the time- 
averaged transition probability (|1.3| ) is likely incorrect far away from MSW resonances. 

An alternative kinetic description based on a production rate obtained from quantum 



field theory has recently been offered [ 44 1 and seems to yield a result very different for 
the rate equation (1.1) in terms of the time-averaged transition probability. However, 
this alternative description focuses on the hadronic contribution near the MSW resonance, 
and as such cannot yet address the issue of the widely separated time scales far away 
from it. A full quantum field theoretical treatment far away from an MSW resonance 
which systematically and consistently treats the two widely different time scales is not yet 
available. 



Thus we conclude that while the result for the rate (|4.2| ) is a direct consequence of 
including the correct transition probability P a ^ s (t) given by ( |3.6[) into the rate equation 
( |1.1| ), our field theoretical analysis of the full neutrino propagator in the medium, and the 
real time evolution of the transition probability, extracted from the full density matrix 



leads us to challenge the validity of the simple rate equation (1.1) with (1.2) to describe 



sterile neutrino production in the early Universe away from an MSW resonance. 
5. Conclusions: 

Motivated by the cosmological importance of sterile neutrinos, we reconsider an important 
aspect of the kinetics of sterile neutrino production via active-sterile oscillations at high 
temperature: quantum Zeno suppression of the sterile neutrino production rate. 

Within an often used kinetic approach to sterile neutrino production, the production 
rate involves two ingredients: the active neutrino scattering rate r aa and a time averaged 



active-sterile transition probability (P a ^ s )[39, 4C, 41, 42, ||] in the case of one sterile and 
one active neutrino. 

Unlike the usual treatment in terms of a truncated 2x2 density matrix for flavor degrees 
of freedom, we study the dynamics of active-sterile transitions directly from the full real 
time evolution of the quantum field density matrix. Active-sterile transitions are studied as 
an initial value problem wherein the main ingredient is the full neutrino propagator in the 
medium, obtained directly from the quantum density matrix and includes the self-energy up 
to 0{G 2 F ). The correct dispersion relations and damping rates of the quasiparticles modes 
are obtained from the neutrino propagator in the medium. We introduce a generalization 
of the active-sterile transition probability from the expectation values of the neutrino field 
operators in the density matrix. 

There are three main results from our study: 

• I): The damping rates of the two different propagating modes in the medium are 
given by 

ri(fc) = r aa (k) cos 2 e m (k) ; r 2 (k) = r aa (k) sin 2 e m (k) (5.1) 

where T aa (k) oc G 2 F k T 4 is the active neutrino scattering rate and 6 m (k) is the mixing 
angle in the medium. The dispersion relations are the usual ones with the index of 



refraction correction [47], plus perturbatively small two-loop corrections of 0(G 2 F ). 
We give a simple physical explanation for this result: for very high temperature 
when 9 m ~ tt/2, v a ~ v% v s ~ v \ an d ~ r aa ;ri <C T aa . In the opposite limit of 



very low temperature and small vacuum mixing angle 9 m ~ , v a ~ v\ ; v s ~ v<i and 
Ti ~ T aa ;r2 <C T aa . Thus in either case the sterile neutrino is much more weakly 
coupled to the plasma than the active one. 

II): We study the active-sterile transition probability P a -> s (t) directly in real time 
from the time evolution of expectation values of the neutrino field operators in the 
density matrix. The result is given by 



-r 1 (k)t + e -r 2 (k)t _ 2e -i(Ti(fc)+r a (fc))t cos ( A £(fc)£) 



(5.2) 

The real time analysis shows that even when T(k) ^> AE(k), which in the literature 
p3| , H is taken to indicate quantum Zeno suppression, the transition probability is 
substantial on time scales much longer than T _1 if T± and T2 are widely separated. 
While the oscillatory interference term is suppressed by the decoherence time scale 



1/T(k) = 2/T aa (k) in agreement with the results of [p4| , 41 1, at very high or low 
temperature this is not the relevant time scale for the suppression of the transition 
probability, which is given by 1/T sm with T sm the smaller between Ti, IV We obtain 
the complete conditions for quantum Zeno suppression: i) 2AE{k)/Y aa <C 1 where 
AE{k) is the oscillation frequency in the medium, and ii) T\ ~ IV This latter condi- 
tion is only achieved near an MSW resonance. Furthermore we studied consistently 
up to second order in standard model weak interactions, in which temperature regime 
the quantum Zeno condition T aa {k) ^> AE(k) is fulfilled. We find that for m s ~ keV 
and 1(T 5 < sin 2d < Kr 3 @, §, g§ the opposite condition, r aa (A;) « AE(k) is ful- 
filled in the high temperature limit T 3> Tmsw ~ 215 MeV, as well as in the low 
temperature regime 3 MeV < T -C Tmsw ■ We therefore conclude that the quantum 
Zeno conditions are may only be fulfilled near an MSW resonance for T ~ Tmsw 
but a firmer assessment of this possibility requires to include the 0{G 2 F ) corrections 
to the index of refraction. 



Ill): Inserting the result ( |5.2|) into the expressions for the time averaged transition 
probability (|1.3|) and the sterile neutrino production rate (|1.2|) yields an expression 
for this rate that is enhanced at very high or low temperature given by equation ( |4.2| ) 
instead of the result ( |1.5| ) often used in the literature. The surprising enhancement 
at high or low temperature implied by (|4.2| ) originates in two distinct aspects: i) the 
assumption of the validity of the rate kinetic equation in terms of a time-averaged 
transition probability with an averaging time scale determined by the decoherence 
scale 2/r aa , and ii) inserting the result ( |5.2j ) into the definition of the time-averaged 
transition probability. The enhancement is a distinct result of the fact that at very 
high or low temperatures the decoherence time scale is not the relevant scale for 
suppression of P a —, s but either 1/Ti or I/T2 whichever is longer in the appropriate 
temperature regime. Our analysis shows that far away from the region of MSW 
resonance, the transition probability reaches its maximum on time scale 1/T(k), 
remains near this value during a long time scale l/r sm ^ 1/r. We have also argued 



that in the early Universe the definition of a time averaged transition probability is 
ambiguous far away from MSW resonances. Our analysis leads us to conclude that 
the simple rate equation ( |1.1| ) in terms of the production rate ( |1,2[ ), ( |1.3|) is likely 
incorrect far away from MSW resonances. 



We emphasize and clarify an important distinction between the results summarized 
above. Whereas I and II are solidly based on a consistent and systematic quantum field 
theory calculation of the neutrino propagator, the correct equations of motion for the quasi- 
particle modes in the medium and the time evolution of expectation values of neutrino field 
operators in the quantum density matrix, the results summarized in III are a consequence 
of the assumption on the validity of the kinetic description based on the simple rate 
equation ( |1 . 1| ) with an effective rate (1.2) in terms of the time-averaged transition proba- 
bility Ql.3| ). The enhancement of the sterile production rate arising from this assumption, 
along with the ambiguity in properly defining a time-averaged transition probability in an 
expanding cosmology in the temperature regime far away from a MSW resonance all but 
suggest important caveats in the validity of the kinetic description for sterile neutrino pro- 
duction in terms of a simple rate equation in this regime. Our analysis suggests that a 
deeper understanding of possible quantum Zeno suppression at high temperature requires 
a reassessment of the validity of the perturbative expansion in the standard model or in 
Fermi's effective field theory. Further studies of these issues are in progress. 
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